NON-NEUTRAL GLOBAL SOLUTIONS FOR THE ELECTRON 
EULER-POISSON SYSTEM IN SD. 
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Abstract. We prove that small smooth irrotational but charged perturbations of a constant back- 
ground are global and go back to equilibrium in the 3D electron Euler-Poissson equation. 
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1. Introduction 

1.1. Presentation of the equation. Consider a three-dimensional plasma composed of fixed 
ions with uniform density, and a gas of moving electrons. This situation can be described by the 
Euler-Poisson equation, which couples a compressible gas to an electrostatic field. Letting u be the 
velocity of the electron gas and p its density, it reads after a simple rescaling 

r dtp + V ■{pu) = 

(1.1) (EP) i p{dtu + u-Vu) = -pVp + pV^ 

( A$ = p-1. 

We took for simplicity the pressure law p{p) = ip^; the analysis is similar for other pressure laws. 
We make the standing assumption that the fluid is irrotational 

(1.2) Vxn = 0; 

this condition is of course conserved by the flow of (|1.1|) . Our aim is to understand the stability of 
the obvious stationary state 

n = , p = I 

under perturbations in p which do not have mean zero (i.e. are not electrically neutral). 

It has been proposed that the non-neutral assumption could have important consequences for 
the asymptotic dynamics of the perturbation. This is made plausible by the following remark: 
consider for each time a ball centered at the center of mass and containing one half of the total 
electric charge. Far away from the ball, the action of the electric field generated by the electrons 
inside this ball is similar to the action of a single point charge at the center of mass of the same 
total electric charge. Recall now that particles moving in a Coulombian electric field experience 
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asymptotically a logarithmic correction to their free trajectory |9]. Based on this analogy, one 



could expect a correction to the linear scattering of the perturbation, which had been shown for 
electrically neutral perturbations (for which the electric field at infinity would be that of a dipole 
and hence integrable along the trajectories). 

It turns out that this picture is not accurate, at least for irrotational perturbations. To the best 
of our understanding, this is due to a combination of the following three facts 

(1) The analogy with point particle is not completely accurate. 

(2) The generated perturbation of the electric field oscillates in time along with the particles. 
However, the nonlinearities come from inertial terms (convection and pressure) and therefore 
oscillate "out of phase" with the density and velocity fields (in other words, the nonlinear 
terms oscillate in a non resonant way with respect to the linear terms). This allows to use 
a normal form transformation and partially cancel their long-time influence (at least cancel 
the influence of the "short-range" part of the nonlinearity involving /?). 

(3) The electric field away from the origin has constant amplitude and oscillates (when there 
is a motion). Its main effect is to periodically repel and attract the particles to the origin. 
However, since its strength is not uniform, the net effect is to create a repulsive force. To 
understand the effect of this force, assume that the center of mass remains at the origin. 
When moving away from the origin, the particles experience a weaker force which takes more 
time to counteract their inertia and invert their velocity, so they move a long distance away 
from the origin; in contrast, when the particles move closer to the origin, they encounter 
a stronger force which sends them back faster and they move a shorter distance towards 
the origin. However, this dispersion effect depends not on the amplitude of the electric 
field (which decays like l/|x|), but on the gradient of this field (which decays like 

and is thus integrable). Therefore, it has no long term effect either and we recover linear 
scattering. 

1.2. Main result. Our main result is that a constant equilibrium of charged electrons is stable, 
even under non neutral (but still irrotational) perturbations. We denote 



JR3 

(which can thus be taken non zero) for the charge of the perturbation. This extends the work of 
Guo [4J who assumed that Q = 0. 

Theorem 1.1. There exists 6q such that if {po,vq) G satisfy 



curl[vo] = 0, II (po - 1,vo)\\hio + \\{po - 1,wo)||h/5,i < 6o, 
then there exists a unique global solution of (jl.ip which converges to equilibrium in the sense that 

\\{p{t)-l,v{t))\\L^ ^Q, ast^+oo. 



Again, as in Guo note the contrast with the result in the absence of electric field in [TT]. This 
theorem suggests that the neutral assumption made in [T, "8] might be removed. Recall that the 
neutral assumption is not necessary either to get small-data/global existence for the Euler-Poisson 
equation for the ions 0, which corresponds to the large-time behavior of the system, but the ion 
case is more transparent since the non-neutral assumption has no implication on the decay of free 
solutions. 

Our analysis relies on works on quasilinear dispersive equations, especially on normal form trans- 
form methods, starting from the work of Shatah [lOj and following recent developments of the 
space-time resonance method in [21 [3l [H [71 [8] . 



(1.3) 




Furthermore, it scatters in a sense which will be made precise in Corollary \4-^ 
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The main consequence of the non neutral assumption Q 7^ is that the solution to the linearized 
equation is no longer integrable in time. However, we remark that its derivative still is. Since the 
quadratic nonlinearities always involve at least one derivative, the main point is then to systemati- 
cally track down the extra decay provided by this derivative term, thus giving a fairly simple proof 
of Theorem II. 1[ 

2. Notations 

We adopt the following notations 

• A < B if A < CB for some implicit constant C The value of C may change from line to 
hue. We note A B ii A < B < A. 

• If / is a function over then its Fourier transform, denoted /, or J-{f), is given by 

m = ^fio = I e-^m dx thus m = ^ J e-q{o d^. 

(in the text, we systematically drop the constants such as ^2j)d/2 since they are not relevant). 

• The Fourier multiplier with symbol m{^) is defined by 

m(iV)/ = [mTf] . 

• The Littlewood-Paley projector P<n, P>n/2 ^iid Pn are defined as the Fourier multipliers 
of symbols 

x(C/iV), (1 - xiUN)) and x(C/(2A^)) - xiC/N) 

for X S C^(M^) such that x(^) = 1 when |x| < 1 and xi^) = when |x| > 2. In what 
follows, sums over capital letters are understood to be over dyadic numbers. 

• The bilinear Fourier multiplier with symbol m is given by 

(2.1) T^[f,g]ix)=^ J e*^«+^)/(0?(r?)m(e,r?)dCdr? = J-^y m(C - r/, r?)/(e - r?)ff(7y) dry. 
We also define Tm to denote an operator "of the form" in the sense that 

Tm [f, g]^ {Tm [/, g] , Tm [7, ff] , [/, g] , Tm [7, } • 



• The Japanese bracket (•) stands for (x) = VT+l?. 

• The Riesz transform is defined as the real operator Rj = \V\~^dj. 

• The Besov spaces are defined by their norms as follows 



E {NrwPNfWh- 



3. Preliminary steps 

In order to investigate the stability of n = 0, /> = 1, using ()1.2p . we introduce the new unknown 
function 

(3.1) a = (V)|Vr^(p - 1) + i|VrMiv[n]. 
The original unknowns can be recovered by the formulas 

p-l = |V|(V)"^9^e[a] and Uj = -Rpm[a]. 
The system (jl.ip becomes 

(3.2) (dt - i(V))a = -lYl ^i^^) [jvl^" + ■ ^^'(^ - ")] - ^ 2^ 1^1 [^J-^" - ■ - ")] • 

j=i ^ ' j=i 



4 PIERRE GERMAIN, NADER MASMOUDI, AND BENOIT PAUSADER 

The above right-hand side is a sum of quadratic terms in a and a: 

RHS = F{a,a) where F{f,g) = Qi{f,g) + Q2{f,g) + Q3{f,g). 
The bihnear operators Qi, Q2, and Q3 are pseudo- products as in (j2.ip whose symbols are hnear 



combinations of the fohowing multiphers 

c\ If I (6 + 6) [6 • (6 + 6)] If i/c \ (^1 + ^2)^ 
-.(6,6) = l^^l (ei) 16116 + 61 = l^^l^^^)7a)(6r-^(^-^^) 

'6 -(6 + 6) 6-6 



?n-t(6,6) = 16 



I6l?7it(6,6) 



16116 + 61 I6II6I 

or their symmetric mp{^2,£,i) and mt{6.2,Ci)- In (|3.3p . mt corresponds to the second term in (|3.2p 
after using that |V| = — Rjdj- 

We now isolate the effect of the electric field at infinity. Actually, for simplicity in our analysis, 
we need to replace Q by Q in p.6p and introduce as follow^ 

/3(t) =^ a{t) - e^*<^>x^ where ^= ^<i^He [a(0)] . 

It solves the system 

= 1 + 11 + III. 
which is forced by e^^^^-'x*^ satisfying (due to (jS.ip ): 

v^'^h'^WB^^ <Q{'^+ tr^^'-'^\ 2<p<3 

(3.5) ||Ve^*<^>X^llBO^ < Q(l + t)"^(^"i\ 2<p<+oo 
uniformly in t > 0, where 

(3.6) g:= ||P<ip(0)||ii 



is a substitute for Q in ()1.3p . 

Fix cj > 2 and > o" + 7. We define our main norm for the global control: 



(3.7) 



Xt ■■= sup 

0<t<T L 



U + t)'^/^||/3(t)||B?,,+ 



Using (|5.ip . we see that, for all T 

(3.8) \\e''^^^f\\xr<\\J\\Y 
uniformly in T. Local existence in Xt follows from energy estimates, therefore we see that Theorem 



11.11 will be a consequence of the following Proposition. 

Proposition 3.1. There exists 6 > such that if (3 G C{[0,T] : H^) satisfies on [0,T] and 

Q+ sup ||/3(i)||H5 <<5, 

0<t<T 

then 

< 



Xt ^ IIPlu;||y + (Q + WWxt) 



uniformly in T. 



^hut one should essentially think of as Q[V| for X ^ nice bump function such that J x^-x = 1. 
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In the proof, we have decided to use Besov spaces instead of more classical spaces in order to 
have a simple access to the estimates in Lemma l5.1[ More elaborate harmonic analysis techniques 
probably allow to replace Besov spaces with Lebesgue spaces. In any case, the difference between 
the two should be thought of as inessential to the proof. Also, we have made some effort to quantify 
the number of derivatives needed and to keep it reasonably low (around 10), notably through a 
"tame" estimat^ on bilinear operators in Lemma 15.11 A slightly more efficient analysis could 
somewhat reduce this number, but to make it close to the number of derivatives in the physical 
energy would require significantly stronger results. 

4. Proof of Proposition 13. II 
We introduce the linear profile 

b{t) = e"'^^^l3{t). 

This natural unknown is only affected by the nonlinearity 

(4.1) dtb = e-**<^> [I + 11 + III] . 
Using Duhamel formula, we see that 

(4.2) /3(t) = e^*<^M/3(0)+AA] 
where A/" is a finite sum of operators of the form 



(4.3) 



Jo 

= J^-^ [ [ e'"^'^^-'^^'^^m{i - ri,ri)C^i{s,i - r])C^2{s,7^)dr]ds 
Jo Jr^ 



for some m{S,i,^2) as in (13.30 . some e = (£1,62) £ ci{t,x),C2{t,x) G {b{t,x),x'^ (x)}, where 

= Id and C~ denotes the complex conjugation and 

(4.4) MCi,C2) = -(6 + 6) + ei(6> + e2{^2)- 

We will obtain bounds which are uniform in e and m. 

In all that follows, the worst term to keep in mind is a variation of e**^^^^'^ • |V|P>ie**<^>6(t) 
where b has to provide both decay and regularity. 

We start with a simple estimate. 
Lemma 4.1. For any choice of c ^ {&)X*^}; ^6 have 

forO<e< 1/100 and 

\\e^'(^^dMt)\\^.-^ < (1 + t)-'/' [Q + mix,]' . 

Proof of Lemma \4-l\ The case c = is trivial. Hence it suffices to treat the case c = b. From 
(j3.3p . (j4.ip . it suffices to show that 

||r^[|V|e^*<^>ci, (V)e^*<^>C2]||^.+3 < (1 + t)-'/'''-'{Q + ||/3|U,)^ 

||r^[|V|e**<^)ci, (V)e**<^>C2]||^.-i < (1 + t)-^/HQ + Mx^f 

for ci,C2 G i^x'^}; where R denotes a Fourier multiplier coming from a symbol of order 0, 
(which we will omit for simplicity of notation) and for fj, £ {1, }; for which it is easily seen 



^tame in the sense that most of the loss of derivative is on the low frequency term 
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(for instance by Coifman Meyer theory) that Holder bounds apply. We first treat the worst case 
using dM]) and (ISTf]) . 

Independently, 



< 



^min(A^5,Ar-i(l + t)-3/2) 



AT 



<(i+r'/'(Q + ii/3iix,)'. 

Interpolating at order (4/10 + e, 6/10 — e), we see that this term is acceptable. 
It only remains to consider, using ()5.2p 

.ao||/3|U.+s<(l + r'/'ll/5|li 



< 



||r^[|V|/3,(V)/3]||^... ^„^„^^ 
This ends the proof. 



□ 



Corollary 4.2. Under the assumptions of the theorem, the solution scatters in H"^ ; namely there 



exists ctoo in such that 



0, as t ^ +00. 



□ 



Proof. This follows immediately from the integrability in time of 

\mt)\\Hs<ii+tr''/'' 

which in turn follows by interpolation from the two bounds in Lemma |4.1[ 

Now, we are ready to prove Proposition 13. 1[ 

Proof of Proposition \3.1[ We start with the dispersive estimate. We claim that under the hypoth- 
esis of Proposition 13. H there holds that 

\Y + {Q + m\xr?- 



(4.5) 



sup (1 + t)5 

0<t<T 



Using ()3.8p . and (j4.2p . it suffices to treat the nonlinear term M. The effect of the nonlinear terms 
Ie,m can be conveniently reformulated through a normal form transformation, which follows from 
a simple integration by parts. Indeed, 

X.,™[ci,C2](t) = -ie~^*<^>T^/^JCfe^*<^>ci(t),C|e^*<^>C2(t)] +ir^/^Jci(0),C2(0)] 

+ i1^e,m/^, [dsCl,C2]{t) + i1^,m/<t>, [ci, dsC2]{t). 



(4.6) 



The first term in (j4.6p can be treated as follows. First, from (j5.4p . Bernstein and Sobolev 
inequalities, we observe that 

E \\P>M/8Tm/^_. [PMe^*<^>Ci(t), P>M/8e^*<^>C2(t)] bf^., 
M>1 



< ^ ||e^*<^)PMCi(t)||^i.oo||e^*<^)c2(t)|| 



M>1 



r>a + 5 
-°10,2 



< (1 + tr'/'{Q + m\x,) M^Mi-'^iic2ii^.+3i 



/5 



M>1 



<(i+r'/'(Q+ii/3iix,) 
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Independently, using (j5.4p with p = oo, we compute that 

M<8 
A/<8 



< 



/2l 



A/<8 



(Q + ll/?llx.)||e^*<^>|V|c2(t)b,V: 



<(l + t)-V2(Q + ||/3||^^)||e^*(v>|v 



C2(t B 



10,2 



If C2 = x*^) this is acceptable; otherwise, using complex interpolation, 



_5 

i; 

and we see that this term is acceptable. Finally, using ()5.3p . 



< (1 + + - < (1 + + m\xr), 



E ll^<A//8^'m/^, [PMe^^^^'^CilO, P>M/8e'*^''^C2(t)] 
A/ 

< 5^ ivf(iV)-+5 ^ ||PMe^*<^>c.(t)||iio||Poe^*<^>|V|c,(t)||z.io 

N<M^O {a,b}={l,2} 

<(i+r'/'(Q+ii/3iix.)' 

For the second term, using (|5.ip . we see that, for k = 1/1000, 

l|e'*<^>T;^M.[ci(0),C2(0)]|b,v, < (l + 0~'/'rrnM.[ci(0),C2 

< (l + t)-6/^|||V|-"ci(0)||^^.+5+2«,20/9|||Vr"C2(0)||^^.+5+2.,20/9 

<(l + t)-^/^{Q+\\b\\x)'. 

We treat the third term in (|4.6p in a similar way, using (jS.ip . Lemma [5. II and Lemma |4. II First, 
we have for the worst term 



Y \\P>M/se'^''^^Ie,m/<p, [P>M/8dsCl,PMC2] Wb-^ ,, 
M 

< V Al + i-5)"'/'l|e''^^^5.ci||^.+3||PA/e'^^^>C2|| ,.ds 
^ Jo " vy,'^ 

< (\l + t - s)-6/5(l + Vmin(M^(l + M)^"^, (1 + s)-^'^'^)ds{Q + 
•^0 M 

<(l + t)-6/5(Q + ||/3|U,)' 
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while the other terms are easier; we now turn to them. Let e be as in Lemma [4.1l and 2^ and (5/2)£ 
be such that 

lie 1 2 e 



then 



2, 2 3' (5/2)e 5 3' 

\\P>M/8e'^^^^Ie,m/<p, [Pm dsCi, P>m/ 8^2] \\b- 



M 



< 



< 



M 



< Al + t - S)-'/' M'\\PMdsC,\\Hl (1 + S] 
M 



-3/10+e 



ds{Q + m\x, 



<(l + t)-6/^(Q + ||/3|U,^ 



and finally, 



M>1 



Y II^Me**^'^^Xe,m/<^JP>4M9sCl,P>4A/C2]||B-Q2 
1 

Y /V + t - s)-'/'ll^'A/r^/^jP>4Me'^<^>a,ci,P>4A/e''^<^>C2] 

V/>1 



IB 



ds 



10/9,2 



< V M'^+^iV r(i + t_s)-6/5||p^a,ci|U2||Poe*'<^>C2|| 5ds 



l<Af<Af~0 



while 



<(i+t)-''/^(Q + ii/3iu,r 

^ ||-PAfe**^'^^Xe^^/<^JP>4A/5sCi, P>4A/C2]||iJ-Q 2 
Af<l 

< [\l + t- s)-^/^||PAfT„/^jP>4A^e^^<''>9,ci,P>4A/e^^<^>C2]||^io/9ds 

M<1 

< 5] M^AO f\i + t-s)-^/' ^||P^5,ci||i2||P^C2|U2 

A//<1 

<(l+t)-6/5(Q+ 11^11^^)2. 



Now, to finish the proof, we appeal to the following result, which follows from a straightforward 
adaptation of the computations in [8]. There exists an energy 

EN.= y [ {\d^p-l)\^ + p\d^v\^ + \\V\-^d''ip-l)\''}dx 



h\<N 



such that, if 



then 



sup ||a(t)||^5 < 1 

0<t<T 



EN{t) ~ ||a(t)||^^ 
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and for any < t < T, 



EN{t)<EN{0)+ \\a{s)\\z'EN{s)ds, 
Jo 



z' < sup(M4 +M3)||PMa||L-. 



M 



Using that 

< min(M^ 1, M"^(M)tt~i)g + min(M~^+^ 1, {Myt~r>)\\(3\\x 
is integrable, we obtain an a priori bound on E^. Tliis ends tlie proof. 

5. Appendix: analytical tools 



□ 



5.1. Linear decay. The standard dispersive estimates for the linear Klein-Gordon equation follow 
from a straightforward application of the stationary phase estimates 



(5.1) 



3 3 3 3 

< tP~2 II fll , . < 

p ,2 



if 2 < p < oo. 



We also use the simple product formula, valid for all 7 G N and smooth functions a, b 

(5-2) \\ab\\H2 ^7 ll«llHjll^llvi/i'" + 

5.2. Boundedness of pseudo-products. We need the following Lemma about boundedness of 
some pseudo-products as defined in (|2.ip . The main feature here is that the extra loss of derivative 
coming from the fact that the phase (ps can be small only gives extra powers of the low frequency 
L. 

Lemma 5.1. Let m be a symbol as in (13. 3p and e G {(it, it)}. Then for any 1 < r < 00 and 

1 < P,Q ^ +00 with 

1 _ 1 1 

r p g ' 

we have Holder's inequality 

(5.3) \\PMf^/^APNa,Pob]\\Lr < H{1 + Lf\\a\\LP\\b\\Li, 

where H = max(M, N, O) and L = min(M, N, O). In particular, we find that for k > 0, 



.+3+.,. + |||V| '^611^1 



"9,2 



Furthermore, for <t > 0, 2 < p < 00, 2 < q < 00, 1 < r < 00, we get 

(5-4) ll-f>M/8^m/0,[^A/a,-P>A//8^]||B?,2 ~ II I ^l^lliJ" 

where < 6* < 7, 7 = 5 i/ M > 1 and 7 = i/ M < 1. 
Proof. Multiplying by a test function, we define 

Im,n,o = {PAic,Tm/<f>^[PNa, Fob]) = jjj^^KM,N,oiyi,y2,y3)a{y2)Ky3)c{yi)dyidy2dy3 



KM,N,o{yi,y2,y3) 



■M' 



Changing variables (^,7?) — ?• r]), we may assume that min(M, 0)/2 < N < 2max(M, O). 
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We compute that 



\Im,n,o\ 



^ KM,N,o{yi,y2,y3)a{y2)b{y3)c{yi)dyidy2dy3, 



< 



E 



S/i-S/2|~A 



|a(y2)||c(yi)|A^6(y2) 



sup B^\KM,N,o{yi,y2,y3) 
\y2-y3\~B 



dyidy2 



< 



"^{ABf sup \KM,N,o{yi,y2,y3) 

\yl-y2\~A, \y2-yshB 



A,B 



\c{yi)\ - Mia- Mb] {yi)dyi 



where M denotes the Maximal functiorH. Using Holder's inequality and the boundedness of the 
Maximal function, in order to prove (j5.3p . it suffices to show that 

(5.5) 



A,B 

We observe that 



Y^{ABfH Hl + L) ^CA,B<1, CA,B= sup I i^M,7V,o (y 1, ^2,^3) I- 

\y1-y2\~A, \y2-y3\~B 



(5.6) 



r.ar.0 mM,N,o{C-V,V) 



< 



HL-^ [e + L 



if L < f 



_n-(|a| + |/3|+2) 



where 9 = |Z(^,r/)|. This follows from the fact that the left-hand side above can be written as a 
linear combinations of terms like 



dl df,^ mM,N,o df <Ae df d'^'^ c\>. 



7i, (^i > 0, 71 H h 7fc = a, (5i H h 5a; = /3 



which is easily seen by induction on |a| + and from the bounds 



< 



•|^|-H|^|H/3| 
|^|-l"l|^|-l/3| [^+1^1 



if L < 1 
if L > 1. 



These bounds follow from elementary but lengthy computations which we omilQ. 
From (I5.6p . we compute that 



-iii-[yi-y2]^-iv [y2 -vsIqc qH '^m,n,o{C - r],r]) 



< 



A^B^\cab\ < sup 

\a\=a; |/3|=fe 

yHLa+b-lj^3-aQ3-b if _L > 1, 

from which (|5.5p follows easily, choosing a, 6 = 2 or 4, 

Y,{ABfH~\Lr^CA,B < Y.{AM/{L)f-\BO/{L)f-' < 1. 

A,B 



-dS^drj 



A,B 



This finishes the proof of ()5.3p . Estimate (|5.4p then follows directly by summing ()5.3p when 
M ~ O > iV. □ 



If r = 1, we have to put the mEiximal function on c instead of on a ■ Mh at the last line above, but this makes 
absolutely no change. 

^Remark however that one can easily obtain worst bounds which would allow for a similar Holder estimate on the 
pseudo-product but with a worst loss in derivative (but still loosing less than 20 derivatives). This would be sufficient 
to obtain the main result assuming more derivatives on the initial data. 
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